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For Anderson localization models, there exists an exact real-space renormalization procedure at 
fixed energy which preserves the Green functions of the remaining sites [H. Aoki, J. Phys. C13, 3369 
(1980)]. Using this procedure for the Anderson tight-binding model in dimensions d — 2,3, we study 
numerically the statistical properties of the renormalized on-site energies e and of the renormalized 
hoppings V as a, function of the linear size L. We find that the renormalized on-site energies e 
remain finite in the localized phase in d = 2, 3 and at criticality {d = 3), with a finite density at 
e = and a power-law decay at large |e|. For the renormalized hoppings in the localized phase, 
we find: In Vl — —-^ — h L'^u, where ^;oc is the localization length and u a random variable of order 
one. The exponent uj is the droplet exponent characterizing the strong disorder phase of the directed 
polymer in a random medium of dimension (d — 1), with t<;(d = 2) = 1/3 and (jj(d = 3) ~ 0.24. At 
criticality (d = 3), the statistics of renormalized hoppings V is multifractal, in direct correspondence 
with the multifractality of individual eigenstates and of two-point transmissions. In particular, we 
measure ptyp ^ 1 for the exponent governing the typical decay In Vl ^ —ptyp InL, in agreement with 
previous numerical measures of atyp = d + Ptyp — 4 for the singularity spectrum /(a) of individual 
eigenfunctions. We also present numerical results concerning critical surface properties. 

I. INTRODUCTION 

In statistical physics, any large-scale universal behavior is expected to come from some underlying renormalization 
('RG') procedure that eliminates all the details of microscopic models. In the presence of quenched disorder, interesting 
universal scaling behaviors usually occur both at phase transitions (as in pure systems) but also in the low-temperature 
disorder-dominated phases. Since the main property of frozen disorder is to break the translational invariance, the most 
natural renormalization procedures that allow to describe spatial heterogeneities are a priori real-space RG procedures 
However, real-space RG such as the Migdal-Kadanoff block renormalizations ^2], contain some approximations for 
most disordered models of interest (these RG procedures become exact only for certain hierarchical lattices [1, Q). In 
this respect, an important exception is provided by Anderson localization Q which has remained a very active field 
of research over the years (see the reviews 0, 0, S, 0, [13, [HI, 1131) • the usual Anderson tight binding model in 
arbitrary dimension d, Aoki [isl . 0, [Tsj has proposed an exact real-space renormalization (RG) procedure at fixed 
energy that preserves the Green functions for the remaining sites (see more details in Section |lT] below). However, the 
numerical results on the RG flows obtained by Aoki thirty years ago were limited to systems of linear sizes i < 16 in 
dimension d = 2 [l^, L < 8 in dimension d — 3 [H, [l3| and to a very small statistics over the samples. The aim of 
the present paper is thus to obtain more detailed numerical results concerning the statistics of renormalized on-site 
energies and renormalized hoppings for Anderson tight-binding model in dimension d = 2, where only the localized 
phase exists, and in dimension d = 3, where there exists an Anderson transition. Our main conclusions are the 
following : (i) in the localized phase in dimension d — 2,3, the statistics of renormalized hoppings is not log- normal 
(in contrast with the conclusions of p^ll4| based on numerics on too small systems), but involves the same universal 
properties as the directed polymer model in dimension 1 + {d — 1), in agreement with [IB, 53, [3 (ii) at criticality, 
' the statistics of renormalized hoppings is multifractal in direct relation with the multifractality of eigenstates (see the 

reviews [13, [ip) and the multifractality of the two-point transmission [l9l. [20ll2]|. 
' The paper is organized as follows. In Section |TT1 we describe the exact renormalization rules for Anderson models 
at fixed energy, and explain the physical meaning of renormalized observables in terms of the Green function. The 
statistical properties of renormalized on-site energies is discussed in Section lHll The statistics of renormalized hoppings 
is studied in the localized phase in Section IIVI and at criticality in Section |Vl Our conclusions are summarized in 
Section rvn 
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II. REAL-SPACE RENORMALIZATION RULES AT FIXED ENERGY 
A. Anderson localization Models 

The renormalization (RG) procedure described below can be applied to any Anderson localization model of the 
generic form 

H = Y,e,\z><z\+Y,V,^j\i><3\ (1) 

i i,j 

where is the on-site energy of site i and where Vij is the hopping between the sites i and j. 

1. Anderson tight binding model in dimension d = 2 and d — 3 

The usual Anderson tight-binding model Q corresponds to the case where 

(a) the sites (z) live on an hypercubic lattice in dimension d 

(b) the hopping Vij is unity if i and j are nearest neighbors (and zero otherwise) 

(c) the on-site energies are independent random variables drawn from the flat distribution 

Pie.) = ^0 < < -) (2) 

The width W thus represents the initial disorder strength. It is known that in dimension d = 1,2, only the localized 
phase exists, whereas in dimension d = 3, there exists an Anderson transition at some critical disorder Wc whose 
numerical value is around (see the review and references therein) 

Wc ~ 16.5 (3) 

2. Power-law Random Banded Matrix (PRBM) model 

The Power-law Random Banded Matrix (PRBM) model is defined as follows : the matrix elements Vij are inde- 
pendent Gaussian variables of zero-mean Vij = and of variance 

= ^ ^ (4) 

where r^.j is the distance between sites i and j . One may consider either a line geometry with rij — \j — i\ or the 
ring geometry of size L (periodic boundary conditions) with 

n \ L 



,i-).3sin(^A_^) (5) 



We refer to our recent works [23|, [2l| for more details and references on the PRBM model. The most important 



property is that the value of the exponent a determines the localization properties 2^1 : for a > 1 states are 
localized with integrable power-law tails, whereas for a < 1 states are delocalized. At criticality a = 1, states become 
multifractal S IM [li, llll 



B. RG rules upon the elimination of one site 

We now consider the Schrodinger equation at a given energy E for an Hamiltonian of the form of Eq. [TJ To 
eliminate a site iq, we use may the Schrodinger equation projected on this site 



j 



(6) 
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to make the substitution 

in all other remaining equations. Then from the point of view of other sites, any factor of the form Vi^ig'4'{io) has to 
be replaced by 

i.e. the hoppings between two neighbors of {q are rcnormalized according to 

rrnew -ir , ^,io^oj' /n\ 

and the on-site energy of each neighbor i of is rcnormalized according to 

e^'^™ = €, + Yl:l°}j:2A (10) 

E - 

These renormalizations equations are exact since they are based on elimination of the variable ipiio) in the Schrodinger 
Equation. The RG rules of Eqs and [TUl have been introduced by Aoki [H, [ij by considering the equations satisfied 
by the Green function. Here we have chosen to derive them in the most elementary way by direct substitution in the 
Schrodinger equation to make obvious their origin and their exactness. 

As stressed by Aoki the RG rules of Egs P and [TUl preserve the Green function for the remaining sites. This 

means for instance that if external leads are attached to all surviving sites, the scattering properties will be exactly 
determined using the renormalized parameters. To get a better intuition of the physical meaning of the rcnormalized 
parameters, it is thus interesting to consider the simplest cases where the disordered system is coupled to only one or 
two external wires as we now describe. 

C. Physical meaning of the renormalized on-site energies 

If one uses the RG rules of Eqs [91 and 1101 until there remains only a single site called A, the only remaining parameter 
is the renormalized on-site energy €a{E). If an external wire is attached to this site A, the scattering eigenstate jV' > 
satisfies the Schrodinger equation 

H\^>^E\i:> (11) 

inside the disorder sample and in the perfect wire characterized by no on-site energy and by hopping unity between 
nearest neighbors. Within the wire, one has thus the plane-wave form 

< xa) = e^'^^^-^-^) -t- re-*'''(^-^-^) (12) 

where the energy E is related to the wave vector k by 

E^2cosk (13) 

The reflexion coefficient r of Eq. [T^j is determined by the ratio 

WyxA) l + r 

that is imposed by the Schrodinger Eq. [TT] projected onto site A. This can be computed in two ways as we now 
discuss. 

1. Solution in terms of the renormalized on-site energy 

In terms of the renormalized on-site energy e^(£'), the Schrodinger Eq. [TT] projected onto site A simply reads 

EiI;{xa) = eA{E)ij{xA) + ij{xA - 1) (15) 
i.e. the ratio of Eq. [14] is directly related to the renormalized on-site energy eA 

R = E-eA{E) (16) 
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2,. Solution in terms of the spectrum of the closed system 

We denote by {En, (fin) the spectrum of the disordered closed system, so that the HamUtonian inside the disordered 
sample reads 

Hsystem ='^En\(f>n X <l>n\ (17) 
n 

In the presence of the wire, the scattering state {i/j > of Eq. [TT] which takes the form of Eq. [T^lin the wire, can be 
decomposed within the disordered system on the (0„) basis 

|V'>=^a„|0„> (18) 



Projecting the Schrodinger Eq. [TT]on < (j)m\ yields the coefficients 

(l>m{xA)ll'{xA - 1) 



E~ En 

In particular at the contact point A, one obtains 



(19) 



i^{xA) ^^an(j>n{xA) = ^^{xa ^ 1) X! ^t^^l]^ (20) 



E — En 

■ft It 

SO that the ratio R of Eq. [14] reads 

\(t)nixAy^ 



R^l^ F- F = Ge{xa,xa) (21) 



R ^ E-E, 

n 

in terms of the Green function of the closed system. 

3. Relation between the on-site energy and the Green function 
In conclusion, the comparison of Eqs [TBI and [211 yields 



1 ^ / N \(t>n{xA)\^ 

- Ge[xa,xa) = 2^ (22) 



E-eA{E) ^' ^' ^ E-En 

i.e. the on site-energy eyi(i?) of the remaining site A is directly related to the Green function Ge^xa, x a) at coinciding 
points. 

D. Physical meaning of the renormalized hoppings 

If one uses the RG rules of Eqs and [TUl until there remains only two sites called A and i3, the only remaining 
parameters the two renormalized on-site energies tAiJ^)^ ^b{E) and the renormalized hoppings VAsiE). 

1. Solution in terms of the renormalized parameters 

In terms of the renormalized parameters, the Schrodinger Eq. [11] projected onto sites A and B simply reads 

E^[xa) = eA{E)i^{xA) + ^{xA - 1) + Vab{E)^xb) 

E^Pixs) = eB{E)^{xB) + ^{xB + 1) + VBA{E)i,{xA) (23) 

If two external wires are attached to A and B the scattering eigenstate \ip > satisfies the Schrodinger Eq. [TT] inside 
the disorder sample and in the perfect wires, characterized by no on-site energy and by hopping unity between nearest 
neighbors, where one requires the plane- wave forms 

'^P^nix < xa) = e''^^^-^-^) + re-^'=(^-^-*) (24) 
i^out{x>XB) =te''^(™) (25) 
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These boundary conditions define the reflection amplitude r of the incoming wire and the transmission amplitude t 
of the outgoing wire. The boundary conditions of Eq. [25] determine the following ratio on the outgoing wire 



i^jxB + 1) 



(26) 



The following ratio 



R 



1p{xA - 1) 



(27) 



concerning the incoming wire can be then computed in terms of the three real renormalized parameters from Eq. 



R^E-€A 

The reflexion coefficient r of Eq. [^5] is then obtained as 

R 



1/2 



-ik 



R 



yielding the Landauer transmission 



T= W = 1 - \r? 



(28) 



(29) 



(30) 



To simplify the discussion, we will focus in this paper on the case of zero-energy E = Q (wave-vector k = tt 12) that 
corresponds to the center of the band. The Landauer transmission then reads in terms of the renormalized parameters 



T{E = 0) = 



[6^(4 + 1) - Vl^eB? + [e| + 1 + Vl^f 



(31) 



For later purposes, it is convenient to rewrite Eqs[23]as a system giving the values 'iP{xa) and 'ijj{xB) at the contact 
points in terms of the values 'ip{xA ~ 1) and 'ijj{x b + ^) the wires as 



iI^{xa) 
iP{xb) 



1 



[E - eA)D 
Vab 



^{XA - 1) + 



Vai 



-MXB + 1) 



iE-eA)iE-eB)D 



with the notation 



D = l- 



{E - eA){E - eB)D 

'^(^^ " 1) + ^ An ^(^^ + ^) 
[E - eB)D 



^ab 



{E-eA){E-eB) 



(32) 
(33) 

(34) 



2. Solution in terms of the spectrum of the closed system 
As above, we denote by [Em 



a,. 



the spectrum of the disordered closed system, (Eq. [T7|) and decompose the 

ts 

(35) 



scattering state on the (0„) basis as in Eq. [THl Projecting the Schrodinger Eq. [TT]on < yields the coefficients 

(I)*^{xa)'4'{xa-1) , 0™(a;s)-0(xB + 1) 



E — E„ 



E — Em 

In particular at the contact points A and B, one obtains 

ip{xA) = ^ an4>n{A) = Ge{xa, xa)iP{xa - 1) + Ge{xb,xa)iI^{xb + 1) 

n 

^^{xb) = ^an(j)n{B) = GE{xA,XB)i^{xA " 1) + G e{x B , X B)i^{x B + 1) 
n 

in terms of the Green function of the closed system 



E ~ E„ 



(36) 



(37) 
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3. Renormalized parameters in terms of the Green function 

In conclusion, the comparison between Eq. [33] and [36] gives the Green functions in terms of the renormaUzed 
parameters 



Gaa 
Gbb 



1 



{E - eA)D 
1 

{E-eB)D 



^-4-8 - jt; 7T§ r^: (38) 

[E - eA)[E - eB)D 



or by inversion the renormahzed parameters in terms of the Green function 

1 



E-eA 

E-€B 



GaaD 
1 

GbbD 



with 



Vab - r^r"" n ^^9) 

'-^AA^-'BB^ 

T/2 

D = l = 1 (40) 

(E-eA)(E-eB) GaaGbb ^ ' 

These relations clarify the physical meaning of the renormalized parameters in terms of the Green functions that are 
usually considered in the literature. 



E. Numerical computations of renormalized parameters 

Periodic B.C. 



i 1 i 1 i 1 i 1 




imp 



FIG. 1: (Color on line) Renormalization procedure in dimension d = 2. The initial state is the tight-binding Anderson model 
on a square lattice of size L^, with periodic boundary conditions in the two directions. Sites are then iteratively eliminated 
using the RG rules of Eqs. |9] and [10] until there remains only the four sites corresponding to the large discs, i.e. there are four 
renormalized on-site energies and four renormalized hoppings at distance L/2 per sample. 

The RG rules of Eqs [1] and [TU] can be followed numerically from the initial condition given by the model of Eq. [1] 
under interest. In the following, we describe the sizes and the statistics over the samples that we have studied for the 
Anderson tight binding model in dimension d = 2 and d = 3 and the PRBM model. 
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FIG. 2: (Color on line) Renormalization procedure in dimension d = 3. The initial state is the tight-binding Anderson model 
on a cubic lattice of size L'^, with periodic boundary conditions in the three directions (here for clarity, the sites of the initial 
model have not be drawn in contrast to Fig. [T] concerning the case d = 2 which is more explicit). Sites are then iteratively 
eliminated using the RG rules of Eqs. [9] and [10] until there remains only the eight sites corresponding to the large discs, i.e. 
there are eight renormalized on-site energies and twelve renormalized hoppings at distance L/2 per sample. 




FIG. 3: (Color on line) Renormalization procedure for the PRBM model with the ring geometry. Sites are iteratively eliminated 
using the the RG rules of Eqs. [5] and [TO] until there remains only the two sites L/2 and L corresponding to the large discs, i.e. 
there are two renormalized on-site energies and one renormalized hoppings per sample. 

1. Anderson tight binding model in dimension d = 2 and d — 3 

For the Anderson tight-binding model described in Section FlI A 11 we have followed numerically the RG rules starting 
from an hypercubic lattice of size L"^ with periodic boundary conditions in all d directions. In each sample, the final 
state that we analyse is an hypercube of linear size L/2, as shown on Figure[l]for d = 2 and on Figure[2]for d = 3 

(i) in dimension d — 2, there are four remaining sites per sample as shown on Fig. (TJ i.e. there are 
four renormalized on-site energies and four renormalized couplings at distance L/2. We have studied the sizes 
L = 12,24,36,48,60,72,84,96,108,120. The corresponding numbers ns{L) of independent samples are of order 
n,{L = 12) = 2.10^ n,{L = 60) = 33.10^, n^(L = 120) = 1150. 

(ii) in d = 3 there are eight remaining sites per sample as shown on Fig. [5J i.e. there are eight renor- 
malized on-site energies and twelve renormalized couplings at distance L/2. We have studied the sizes L = 



8 



4,6,8,10,12,14,16,18,20,22,24,26,28,30. The corresponding numbers ns{L) of independent samples are of order 
n,{L = 4) = 10^ ns{L = 10) = 6.10"*, n,{L = 20) = 400 and n,{L = 30) = 24. 



2. Power-law Random Banded Matrix (PRBM) model 



For the PRBM model described in Section fll A 21 we have followed numerically the RG rules up to the final state 
shown on Fig. [3] containing only the sites L/2 and L, i.e. in each sample, there are two renormalized on-site energies 
and one renormalized coupling. We have studied rings of sizes 50 < L < 1800 with corresponding statistics of 
10.10^ > ns{L) > 2400 independent samples. 



III. STATISTICS OF RENORMALIZED ON-SITE ENERGIES 



A. General properties 



We find that the renormalized on-site energies remain finite in all phases (localized, delocalized, critical), and that 
the histograms 'pL(e) corresponding to various system sizes L converge towards some stationary distribution 'Poo(e) 
that present the following common properties : 

(i) 'Poo(e) is symmetric in e ^ — e (as the initial condition of Eq. [2|) 

(ii) 7^oo(e) has a finite density 7^00(0) at its center e = 0. After the change of variables to y = In |e|, this corresponds 

to 

Poo(2/ = ln|e|) oc |e| = e« (41) 

(iii) For e ±cx), Voo{^) presents the following power-law decay 

^-(e), (42) 



C— >±CXD g 

After the change of variables to ?/ = In |e|, Eq. l42l corresponds to 



P^{y = \n\e\) cx ^ = e"^ (43) 



The origin of the power- law of Eq. 1421 even when one starts from a bounded distribution in e as in the tight-binding 
Anderson model (see Eq. [2]), can be understood from the form the RG rule of Eq. [TOl which reads at zero energy 
E = Q 



Vi,io ^io,i 



(44) 



During the first steps of renormalization where the hoppings V are finite, very large renormalized on-site energies are 
generated when the eliminated on-site energy is very small. The finite density of 7^(6^^) at = yields the power- 
law decay of Eq. [42lvia the change of variable e™™ ~ — 1/eio using the standard formula P„eu,(e"'^"')de"'^"' = V{eig)d€ig. 
In the remaining of this section, we present the histograms we have measured in various cases. 



B. Results for the square lattice in dimension d = 2 

On Fig. m we show the histograms of the logarithm of the absolute value of the renormalized on-site energy e for 
various sizes 12 < L < 120 : apart from the cut-offs in the tails imposed by different statistics over the samples, 
these histograms coincide. This shows that the convergence towards the stationary distribution Poo(e) is quite rapid 
: starting from the initial condition of Eq. [21 our results for the smallest size I = 12 have already 'converged' towards 
the final- and very different- distribution of Figure S) On Fig. 31 the slope of the left tail is of order -1-1 in agreement 
with Eq. [4TI and the slope of the right tail is of order —1 in agreement with Eq. |43l 
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In P 




FIG. 4: (Color on line) Statistics of renormalized on-site energies e for the Anderson model on the square lattice in dimension 
d — 2 for sizes 12 < L < 120 (disorder strength W = 40) : the histograms of In |e| are identical (apart for the cut-off imposed 
by the statistics over the samples) The left and right slopes of value unity corresponds to Eqs 31] and 1431 
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FIG. 5: (Color on line) Statistics of renormalized on-site energies e for the Anderson model on the cubic lattice in dimension 
d = 3 for sizes 4 < L < 20 (a) Histograms of In |e| in the localized phase (W = 40) (b) Histograms of In |e| at criticality 
(Wc = 16.5) 



C. Results for the cubic lattice in dimension d — 3 



Our data for the Anderson tight binding model in d = 3 are shown on Fig. [5] : both in the localized phase and at 
criticality, the convergence in L towards the stationary distribution Vooi^) is still rapid, and the measured tails are 
again in agreement with Eqs 1411 and 1431 It turns out that for a given disorder value, our numerical results concerning 
■Poo(e) seem to coincide for d — 2 and d — 3 (Fig. [Hand Fig. [5] (a) corresponding to Wd — 40) : the reasons of this 
coincidence are not clear to us, since the initial coordinence of sites clearly depends on the dimension d. 
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FIG. 6: (Color on line) Statistics of renormalized on-site energies e in the PRBM model of parameter b = 0.1 (a) The 
histograms of In je| at criticality a = 1 for various sizes L = 100, 200, 400, 600, 800 are identical (apart for the cut-off imposed by 
the statistics over the samples), (b) Comparison of the stationary distributions in the localized phase (a = 1.4), at criticality 
(o = 1) and in the delocalized phase (a = 0.6). The left and right slopes of value unity corresponds to Eqs 1411 and 1431 



D. Results for the PRBM model 



The properties found above for Anderson tight-binding models seem to be valid for more general Anderson models 
of the form of Eq. [TJ As an example, we show on Fig. [6] our data concerning the PRBM model described in section 
III A 21 The histograms of renormalized energies converge rapidly towards their limit. The stationary distribution 
presents the tails of Eqs BD and 1151 in all phases (locaHzed, critical, delocahzed). 

E. Consequences 

In conclusion, the renormalized on-site energies remain finite random variables in all phases (localized, critical, 
delocalized). As a consequence, the behavior of the two-point Landauer transmission of Eq. [31] is determined by the 
properties of the renormalized hoppings 

(i) in the delocalized phase, both the renormalized hopping and the two-point transmission will remain random 
finite variables. 

(ii) in the localized phase and at criticality where the two-point transmission decays with the distance, its decay 
will be directly related to the decay of the renormalized hopping via 

In T{E 0) ~ In Vlg + finite (45) 

In the following, we discuss the statistics of renormalized hoppings in the localized phase and at criticality, in 
relation with the statistics of two-point transmission. 

IV. STATISTICS OF RENORMALIZED HOPPINGS IN THE LOCALIZED PHASE 

A. Universality class of the directed polymer in a random medium 

In dimension d ^ 1, the transfer matrix formulation of the Schrodinger equation yields a log- normal distribution 
for the Landauer transmission [l^, 

lnT('^=^) oc -A + Li/2„ (46) 
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The leading non-random term is extensive in L and involves the locahzation length ^/oc ■ The subleading random term 
is of order L^/^, and the random variable u of order 0(1) is Gaussian distributed as a consequence of the Central 
Limit theorem. Although it has been very often assumed and written that this log-normal distribution per sists in 
the localized phase in dimension d— 2,3, theoretical arguments Tt'I and recent numerical calculations Tl8| are in 
favor of the following scaling form for the logarithm of the transmission 



(47) 



where the exponent Lo{d) depends on the dimension d and coincides with the droplet exponent characterizing the 
strong disorder phase of the directed polymer in a random medium of dimension 1 -I- (d — 1) (see the review [29| on 
directed polymers). The probability distribution of the rescaled variable u is not Gaussian but is determined by the 
directed polymer universality class (see [ld | where its distribution in = 2 is shown to coincide with the exactly 
known Tracy- Widom distribution for the directed polymer in 1-1-1). 

The arguments in favor of the same universality class can be decomposed in two steps [l^ [13, [13 ■ 

(i) in the localized phase of Anderson localization in dimension rf, the transmission decays exponentially with the 
length, and thus directed paths completely dominate asymptotically over non-directed paths. In dimension d = 2, 
the dominance of a narrow channel can be seen on Figs 10 and 11 of Ref [30| . 

(ii) these directed paths of the Anderson model have weights that are random both in magnitude and sign, but it 
turns out that the directed polymer model which is usually defined with random positive weights (Boltzmann weights) 
keeps the same exponents in the presence of complex weights (see section 6.3 of the review [29|). 

In conclusion, from the relation of Eq. 1451 we expect that the renormalized hoppings will present the same statistics 
as the Landauer transmission of Eq. [47l 



InVi 



L 

^loc 



(48) 



To check this relation, we have measured the averaged value and the variance of the logarithm of the renormalized 
hoppings in dimension d = 2, 3. 



B. Results for the square lattice in dimension d — 2 
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FIG. 7: Statistics of the typical renormalized hopping in d = 2 where only the localized phase exists (the data shown correspond 
to the disorder strength W = 40) (a) Typical exponential decay : In Vl is linear in L and the slope represents the inverse of 
the locahzation length £,ioc (Eq. HIJ. (b) The fluctuation term A(ln Vl) grows as (Eq. HH with uj{d = 2) ~ 0.33. 



In dimension d — 2, only the localized phase exists. On Fig. [7] (a), we show the typical exponential decay 
corresponding to a finite localization length £^ioc in Eq. 051 On Fig. [7] (b), we show the amplitude A(lnVL) of the 
random term in Eq. |48] : the three parameters fit A(ln V^) = aoL'^*-^^^-' -I- oi yields the value 



uj{d ■ 



0.33 



(49) 



12 




FIG. 8; (Color on line) Histograms of the logarithm of the renormalized hopping InVi for various lengths L in the localized 
phase (for the disorder strength W = 40 (a) in dimension d — 2 (b) in dimension d = 3 



in agreement with the exact result [3l|, |33, 123, [34 



iODp{l + l) = - 



(50) 



for the directed polymer in a random medium of dimension 1 + 1. On Fig. [5] (a), we show the histograms of (In V^,) 
for various sizes L : as L grows, the maximum moves linearly while the width grows as L'^ . 



C. Results for the cubic lattice in dimension d = 3 
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FIG. 9: Statistics of the typical renormalized hopping in d = 3 in the localized phase {W = 40) (a) In Vl as a function of L 
: the slope represents the inverse of the localization length ^joc (Eq. 148 |l (b) The fluctuation term A(lnyL) grows as (Eq. 
148)) with tj(d = 3) ~ 0.24. 

In dimension d = 3, the locahzed phase exists in the domain W > Wc — 16.5 for the disorder strength (Eq. [3]). 
The data shown on Fig. [9] corresponds to the disorder strength W — 40. The histograms of (InVz,) for various sizes 
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L are shown on Fig. [8] (b). On Fig. [9] (a), the typical exponential decay found corresponds to a finite localization 
length £^ioc in Eq. US] On Fig. [5] (b), the amplitude A(ln V^) of the random term in Eq. [35] can be fitted by the form 
A(lnyL) = aoL"(''=3) + ai that yields the value 

w(d = 3)~0.24 (51) 

in agreement with the measures of the droplet exponent uj npi l + 2) ~ 0.244 obtained for the directed polymer in a 
random medium of dimension 1 + 2 in various Refs [H, IsgI. ISTI. [38l [39l | . 



V. STATISTICS OF RENORMALIZED HOPPINGS AT CRITICALITY 



A. Expected multifractal statistics 



At criticality, the statistics of the two-point transmission is multifractal [l9|, |20|, |2l| : the critical probability 
distribution of the two-point transmission Tl takes the form 

Prob (Tl - L-'') dT cx L^^^Uk (52) 

where the multifractal spectrum $(«) exist only for k > (as a consequence of the physical bound 7l < 1) and is 
related to the singularity spectrum f{a) of eigenfunctions via 

$(k > 0) = 2 /(a = d+ |) -d (53) 

At criticality the decay of the two-point transmission is directly related to the decay of the renormalized hopping 
via Eq. 1451 As a consequence, what is known about the statistics of the two-point transmission at criticality can 
be translated for the renormalized hoppings. The probability distribution of the renormalized hopping Vl at scale L 
takes the form 

VT:oh{\VL\-- L'P)dV (x L^'^PUp (54) 

L—^OQ 

where 

i7(p > 0) = $(2p) = 2[/(a = d + p) -d] (55) 
In particular, the typical exponent ptyp characterizing the typical decay 



InVi ~ -ptyplnL (56) 

L — >+oc 

is related to the typical exponent Ktyp of the two-point transmission and to the typical exponent atyp of the singularity 
spectrum via 

Ptyp = ^ = atyp - d (57) 



B. Results for the cubic lattice in dimension d — 3 

The typical exponent ptyp of Eq. [55] is measured from the data of Fig. [TU] (a) 

Ptyp ^ 1.05 (58) 

Via Eq. 1571 this value is in agreement with the numerical measures of order (40l . |4l| 

atyp ~ 4 (59) 

for the exponent atyp concerning the singularity spectrum of eigenfunctions. 

Of course, beyond this typical exponent, one could in principle extract from our numerical data, results on the whole 
multifractal spectrum. However, our numerical means in d = 3 are limited to rather small sizes and small statistics 
(see the section Hi E II) in comparison with the exact diagonalization calculations of Refs f4l|. As a consequence, our 
numerical results seem sufRcient to measure the correct typical exponent, as shown above, but we believe that they 
are not sufficient to measure correctly the rare events that are necessary to obtain a reliable multifractal spectrum. 
It may be that in the future, more 'professional numericians' will be able to transform the present renormalization 
approach into a competitive numerical method to measure the multifractal spectrum, but this is clearly beyond our 
numerical means. 
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FIG. 10: Statistics of the typical renormalized hopping in d = 3 at criticality {W = 16.5) (a) InVi, as a function of InL for 
bulk sites (see Fig. [5]) : the measured slope ptyjj — 1.05 represents the typical exponent of Eq. 1561 (b) InV^"'-^ as a function 
of InL for surface sites (see Fig. Ill[) : the measured slope pty'p'^ — 1.6 represents the typical exponent of Eq. 1601 

C. Renormalized hopping between two surface points in dimension d = 3 




FIG. 11: (Color on line) Renormalization procedure in dimension d = 3 to measure the renormalized hopping Vl^^ between 
two boundary sites. The initial state is the tight-binding Anderson model on a cubic lattice of size L^, with periodic boundary 
conditions in two directions x and y and free boundary conditions in the third direction z. Sites are then iteratively eliminated 
using the the RG rules of Eqs. |5] and [TO] until there remains only the two surface sites corresponding to the large discs. 

At criticality, points lying on the boundaries are characterized by a specific multifractal spectrum fsurf (a) , different 
from the bulk spectrum f{a) .42 , 43. . 44] • These surface critical properties are particularly interesting in Anderson 
localization models where it is more natural to attach leads to boundary sites rather than bulk sites. We have thus 
considered the renormalization procedure depicted on Fig. [Tl]to measure the statistical properties of the renormalized 
hoppings between two surface points. The typical behavior shown on Fig. [TO] (b) 



(60) 
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corresponds to an exponent of order 

pZ^ ^ 1-6 (61) 

clearly distinct from its bulk analog of Eq. [58l 

Via Eq. [571 we expect that this value corresponds to 

a^/ = d + pZ/ (62) 

for the typical exponent alyp^ of the surface singularity spectrum of eigenfunctions. In contrast with the bulk case, 

we are not aware of any direct measure of cttyp'^ the literature to make some comparison. As explained at the end 
of section rVB[ we believe that our numerical means are not sufficient to measure correctly the rare events to obtain 
the full multifractal spectrum around this typical value. However, we expect that our result for the typical exponent 
is reliable (as shown above for the bulk case), and will be confirmed in the future whenever the surface multifractal 
spectrum will be measured via the powerful exact diagonalization techniques of [4l'| . 

From Eq. 1571 we also expect that the two-point transmission in d = 3 between two surface points involves the 
typical exponent 

'^typ^ — "^Ptyp^ — 3-2 (63) 



D. Results for the PRBM model 



For the PRBM model, we have studied in detail the multifractal properties of the two-point transmission in our 
previous works [13, HH ■ Since the statistics of renormalized hoppings can be directly deduced from them via Eqs [55l 
we refer the interested reader to [lOj HH where we have measured multifractal spectra at criticality a — 1 for various 
values of the parameter h. 



VI. CONCLUSION 



In th is p aper, we have revisited the exact real-space renormalization procedure at fixed energy proposed by Aoki 
[3 [3 foi' Anderson localization models. We have presented detailed numerical results concerning the statistical 
properties of the renormalized on-site energies e and of the renormalized hoppings 1^ as a function of the linear size 
L for the Anderson tight-binding models in dimension d = 2 where only the localized phase exists, and in dimension 
d = 3 where there exists an Anderson localization transition. Our main conclusions are the following : 

(a) the renormalized on-site energies e remain finite in the localized phase in d = 2, 3 and at criticality (d = 3), 
with a finite density at e = and a power-law decay 1/e^ at large |e|. 

(b) in the localized phase in dimension d — 2,3, the statistics of renormalized couplings belon gs t o the universality 
class of the directed polymer in a random medium of dimension 1 -f (d — 1), in agreement with jl6l . [TtI . [isj . 

(c) at criticality, the statistics of renormalized hoppings V is multifractal, in direct correspondence with the mul- 
tifractality of individual eigenstates and of two-point transmissions. In particular, our measure ptyp — 1 for the 
exponent governing the typical decay In Vl — —ptyp InL, is in agreement with previous numerical measures of 
o^typ = d + Ptyp ~ 4 for the singularity spectrum /(a) of individual eigenfunctions. We have also measured the 
corresponding critical surface properties. 

In conclusion, we have shown that the large scale properties of Anderson localization models actually emerge from 
the simple real-space RG rules of Eqs [51 and [TUl which preserve exactly the Green functions of the remaining sites. 
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